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Abstract 



The inclusive S"— wave charmonium production rates in B decays are con- 
sidered using the Bodwin-Braaten-Lepage (BBL) approach, including the 
relativistic corrections and the color-octet mechanism suggested as a pos- 

I / 

r~| , sible solution to the -0 puzzle at the Tevatron. We first consider relativistic 

and radiative corrections to J/ijj e~^e~ and J/tp ^ Light Hadrons (LH), 
^ \ in order to determine two nonperturbative parameters, {J/'ilj\Oi(^Si)\J/'ip), 

' {J / 'ij)\Pii^ Si)\J / ij)) , in the factorization formulae for these decays. Using 

^ ■ these two matrix elements and including the color-octet cci^Si) state con- 

tribution, we get a moderate increase in the decay rates for B decays into 
J/il) (or tp') + X. Our results, B{B J/i/j (or V') + X) = 0.58 (0.23)% 
for Mb = 5.3 GeV, get closer to the recent CLEO data. As a byproduct, we 
prefer a larger decay rate for r]c — > LH compared to the present data. 



Typeset using REVT^ 



* pko@phyb . snu . ac . kr 
^ j ungil@phyy. snu . ac . kr 
^hssong@phyy.snu.ac.kr 



1 



I. INTRODUCTION 



It has been commonly believed that inclusive production rates of a heavy quarkonium 
state in various high energy processes can be adequately described by perturbative QCD 
(PQCD) within the color-singlet model [||]. However, recent observations of inclusive J/tj) 
and ip productions at high pt at the Tevatron suggest that a new mechanism is called for 
beyond the color-singlet model 0. For the J/ip production, the lowest order subprocess 
comes from parton fusions, which are smaller than the data by more than an order of 
magnitude. If one includes the gluon fragmentations into J/tp and into Xcj(l-P) states 
followed by Xcj(l-P) ^ J lip + (which is the next-to-leading order in a^), one gets higher 
theoretical estimates that still underestimates the experimental yield by factor of 3 ~ 5. For 
the ip production, the situation is even worse. Even with the gluon fragmentation included, 
the theoretical production rate falls below the data by factor of ~ 30 or so. 

In order to resolve this puzzle, basically two scenarios have been suggested up to now : (i) 
existence of new charmonium states above the DD threshold, which can decay into J / ip and 
ip with appreciable branching ratios 0- and (ii) importance of the gluon fragmentation 
into a pointlike color-octet S"— wave cc state, and its subsequent evolution into ip 

0. Both scenarios are quite intriguing in a sense that they call for new elements of physics 
within the standard model, new spectroscopy or new production mechanism for charmonium 
states. It would be useful to explore and test these suggestions in places other than pp 
colliders such as the Tevatron. In Ref. |^, one of us has explored the consequences of the 
first scenario, finding that these hypothetical Xcj(2-P) states should be observed at the level 
of (0.3 — 0.5)% (in branching ratio) in the decay channels of i? — >• Xc,j=i,2(2P) + X (followed 
by Xc,j=i,2 —*■'?/'+ 7 with ^ 10% in branching ratio, if the first scenario is to work). 

In this paper, we explore the consequences of the second scenario in the inclusive decays 
of B mesons into a S*— wave charmonium. It is well known that the lowest order results 
in the heavy quark velocity and the strong coupling constant for r(_B J/ip + X)direct 
is smaller than the CLEO data by factor of ~ 3 The situation does not get better 
even if the next-to-leading order corrections in as are included in the nonleptonic effective 
weak hamiltonian for B decays In Ref. [Q, only the color-singlet contribution has been 
included, since the color-octet contributions are higher order in f ^, hence suppressed relative 
to the singlet contributions. However, the Wilson coefficient for the color-singlet contribution 
is much suppressed compared to that for the color-octet contribution. Therefore, the color- 
octet contribution, being suppressed by v^, may be numerically important, because of the 
larger Wilson coefficient. This is similar to the case of the gluon fragmentation into Xcj(l-P) 
states [|10|, for which the color-octet contribution is lower order in compared to the color- 
singlet contribution, whereas both of them are of the same order in v"^. Also, in the case of ip 
production through the gluon fragmentation into a color-octet cc state, the octet contribution 
is higher by compared to the color-singlet contribution, but this is compensated by the 
large short distance factor l/a^ compared to the color-singlet contribution 0. 

In Sec. |T|, the 5— wave charmonium production rates in B decays are calculated in 
the framework of Nonrelativistic QCD (NRQCD) |]1T| including the relativistic corrections 
and a color-octet [cc(^S'i)] contribution which is of 0{v'^) compared to the nonrelativistic 
limit. The results contain three nonperturbative parameters, (0|Of^(^S'i)|0), {Q\Pi Si)\Qi) 
and (0 1 (^5*1)10). Among these three parameters appearing in the heavy quarkonium 
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productions, the first two color-singlet matrix elements can be related another parameters, 
{H\Oi(^Si)\H) and {H\Pi(^Si)\H), which enters in the heavy quarkonium decays in the 
vacuum saturation approximation |11|], via 



(O|0f |0) ^ (2J + 1) {H\0^\H) (l + 0{v')) . (1) 
In Sec. inrai , the two parameters {H\Oi{^Si)\H) and {H\Pi{^Si)\H) (with H = J/tp,^') 



are determined by analyzing the decays of J/ijj and ^p' into light hadrons (LH) and e"^e" 
Implications of this analysis on the decays of 7]^ into light hadrons and 77 are discussed. Our 
result prefers larger decay rate for rjc —>■ LH. In Sec. |III B| , we present numerical estimates 
for S"— wave charmonium production rates in B decays, using the factorization formulae 
obtained in Sec. |l| and (O|0f (3^i)|0), {0\Pl^{^Si)\0) obtained in Sec. [ITTA| . We also discuss 



the polarization of the J /ijj in B decays. In Sec. |V|, the results are summarized, and possible 
improvements of the present work are speculated. 



II. HIGHER ORDER CORRECTIONS 

The effective Hamiltonian for b ccq ( with q = d,s) is written as 



G 



F 



Heff = VcbV*g 



2C — C_ 

— ^ C7m(1 - 75)c qY{l - l5)b 



+ (C+ + C^) C7m(1 - j,)T^c q^il - 75)T"6] , (2) 

where C±'s are the Wilson coefficients at the scale /i ~ Mb. We have neglected penguin 
operators, since their Wilson coefficients are small and thus they are irrelevant to our case. 
To leading order in as{Mb) and to all orders in as{Mb) hi{Mw/Mb), the above Wilson 
coefficients are 

C+{Mb) ^ 0.87, C_(Mfc) ^ 1.34. (3) 

According to the factorization theorem for the S*— wave charmonium productions in B 
decays, one has [Q 

Tib - J/^ + X)= (OlQi'^'^^^OlO) _ (^-)^(35^) + (4) 

nb - + X) = (^^orc^'^m ^^^^ _^ ^ 

in the nonrelativistic limit, where Fi are rates for hard subprocesses of b quark decaying 
into a cc pair with suitable angular momentum and vanishing relative momentum in the 
color-singlet : 

fi(6 ^ (c-ches,) + s,d) = - (^1 + fo, (6) 

f 1(6 ^ {cc)^CSo) + s,d) = (2C+ - fo, (7) 
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with 



^o^\V.,nTjt M^MJl--^ . (8) 



The operator (^^^^Sj) is defined in terms of heavy quark field operators in NRQCD 
Its matrix element {O\0i {^^^^Sj)\0) contains the nonperturbative effects in the heavy 
quarkonium production processes, and is proportional to the probability that a cc in a 
color-singlet S*— wave state fragments into a color-singlet 5— wave cc bound state such as 
a physical J/ip, rjc or ip . It is also related to the matrix element {H\Oi(^^^^ Sj)\H) as in 
Eq. (1), and also with the nonrelativistic quarkonium wavefunction as follows : 

{0\Oi/^CS,m ^ 3 {JMO^CS,)\J/^) ^ (A) \R40)f , (9) 

in the nonrelativistic limit. Similar expressions hold for the case of nonperturbative matrix 
elements appearing in the rjc productions and its decays : 

{o\orcsom ^ {vcm'So)\vc) ^ \r,m\" ■ m 

Note that dependence on the radial quantum numbers n enters through the nonperturbative 
parameters, {0\O^ Si)\0) . 

Using the leptonic decay width of J/ip and ip , one can determine 

{J/^\OiCSi)\J/^) ^ 2.4 X 10-^ GeV^ (11) 
{tlj'\OiCSi)\^') ^9.7 xlQ-^ GeV^ (12) 

in the nonrelativistic limit with as{Mc) = 0.27. Q Also, to the lowest order in v'^, one has 
{Vc\Oi(^So)\r]c) = {J/ iplOii^^ Si)\J / ip) because of heavy quark spin symmetry From 



these expressions with Mf, ^ 5.3 GeV, one can estimate the branching ratios for B decays 
into J/i/j + X and ip' + X : 

B{B ^ J/i; + X) = 0.23%, (0.80 ± 0.08)%, (13) 
B{B ^r], + X) = 0.U%, (<0.9% (90%aL.)), (14) 
B{B ^ ij' + X) = 0.08%. (0.34 ±0.04 ±0.03)%. (15) 



^We follow the notations in Ref. and will not give explicit forms for these dimension-six 
operators in this paper. 

^The radiative corrections in has not been included here for consistency. To be consistent with 



the velocity counting rules in the NRQCD in the Coulomb gauge for the heavy quarkonia |11], one 
has to include the relativistic corrections as well, since v ~ as{Mv) in heavy quarkonium system. 
If one includes the 0{as) radiative corrections to J/ip ^ without relativistic corrections, one 
gets a larger (0|O/^^(^5i)|0) compared to the lowest order result, Eq. (11) : {^\0'l^^ Si) 
4.14 X 10~^ GeV^. Relativistic corrections gives a further enhancement. See Eq. (48) below. 
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The recent data from CLEO are shown in the parentheses, where the cascades from 
B — i> Xcj(l-P) + ^ followed by Xcj J/ip + 1 have been subtracted in the data shown. 
In view of these results, we may conclude there are some important pieces missing in the 
calculations of decay rates for B — > (cc)i(^S'i) + X using the color-singlet model in the 
nonrelativistic limit. 

In Ref. P], it was noticed that the inclusion of color-octet piece, which is often ne- 
glected in the previous studies of the charmonium production in B decays, is mandatory 
in order to factorize the amplitude consistently without any infrared divergence in case of 
B decays into the P— wave charmonium. This also leads to nonvanishing decay rates for 
B {he, XcO) Xc2) + X, all of which vanish in the color-singlet model. 

In view of this, we first estimate the color-octet contributions to S — > J/ip+X, motivated 
by the suggestion that the color-octet mechanism might be the solution to the ip puzzle at 
the Tevatron. Although it is of higher order in (~ 0{v^)), it can be important in the 
case of the inclusive B decays into J/ip + X, since the Wilson coefficient of the color-singlet 
part is suppressed compared to that of the color-octet part by a factor of ~ a^- (In Eq- (2), 
(2C+ - C_) ^ 0.4, and (C+ + C_) ^ 2.20. ) Now, it is straightforward to calculate the 
contribution of {cc)s{^Si) to B ^ J/ip + X : 



TiB icc)sCS,) +X^J/^ + X)= ^"'^^^m/'^'"^ + + ^ 1 ^''^ 



TiB (ccU'So) +X^r^e + X) = ^ ^"'^^^f^f °^ (C^ + C.f f,. (17) 

(Similar expression holds for the B ^' {or rj') + X except that (O|0f (^^^OjO) and MjMb 
should change appropriately in order to account for the phase space effects.) 

Here, a new nonperturbative parameter (0|O8(^>S'i)|0) comes in, which creates a (cc) 
pair in the color-octet state, projects into the subspace of states which contain J/ip in the 
asymptotic future, and then annihilates the (cc) pair at the creation point. The matrix 
element of this operator is proportional to the probability that the color-octet (cc) 8(^5*1) to 
fragment into the physical J/tp state in the long distance scale. This type of a color-octet 
operator was first considered in the gluon fragmentation function into P— wave charmonia in 



Ref. ||10[, and then in the gluon fragmentation into ip to solve the ^jJ puzzle at the Tevatron. 

Braaten and Fleming fixed (0|of (^S^jO) to be 4.2 x 10"^ GeV^ (for M, ^ 1.5 GeV), in order 
to fit the total cross section for the inclusive ip production cross section at the Tevatron, 

and found that this value of (OjOg (^5*1)10) yields the pt spectrum for the ip' production 
which nicely agrees with the measured shape. Then, Cho and Leibovich have performed 
a complete analysis for the color-octet contribution to Upsilon and Psi productions at the 
Tevatron for both low and high pt regions |T3|. Their results are 

(0|O8/^('^i)|0) = 1.2 X 10-2 GeV^ (18) 

{O\0t' CSi)\0) = 7.3 X 10-^ GeVl (19) 
Taking the ratio between the color-octet and the color-singlet contributions, one gets 
T{B^{c-cU'S,)+X-.H + X) ^ 3{0\O^i'Sm jC^ + C^)' 

T{B^{cc),{^Si)+X ^H + X) 2(0|Of(3Si)|0) (2C+ - C_)2 ^ 

= 0.76 (1.14) for H = J/tP {tp'), (21) 
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Thus, we find that the color-octet {cc)s{^Si) contributions to B ^ J/ip (or ip') + X are 
about 76 % (114 %) of the color-singlet contributions in the nonrelativistic hmit. 

There is another subprocess in the lower order in compared to the color-octet contri- 
bution considered above : the relativistic corrections to the color-singlet component which 
is an order of O(f^) compared to (4), (5). Extending the Feynman rule in the presence of a 
heavy quarkonium , we can derive that the relativistic corrections to B ^ J/ip (or rj^) +X 
can be written as the following factorized form : 

Tib - J/^ + X) = _ (0|-P'^"^^i)|0) f ^ (3^^) + X), (22) 
Tib ^Vc + X) = J^^L^pR f,(6 (c-c), CS,) + X). (23) 

Here again, one can use the vacuum saturation approximation, Eq. (1) : 

{0\PfC^+^Sj)\0) ^ (2J + 1) {H\PiC^+^Sj)\H). (24) 

The latter is related with the spin-weighted average of the S*— wave heavy quarkonium wave- 
functions in the following way : 



{nS\P,\nS) = i^ns ^'Rns) ^ ^^^2)\ ^ (25) 

27r ^ ^ 



RnS is the spin- weighted average of the 5— wave wavefunctions |11| : 



RnS = ^ (3-Rv + Rvo) ■ (26) 

In order to estimate the relativistic corrections, we need one more nonperturbative matrix 
element, (0|P"'^(^S'i)|0) or Re {R^s V^-R^s). This is not available in the current literature 
now, and we will determine this parameter as well as (0|Oi('^S'i)|0) in the following section. 
However, we note that the relativistic corrections make the decay rates for the S"— wave 
charmonium productions in B meson decays decrease because (0|P"'^(^S'i)|0) > 0. 



III. NUMERICAL ANALYSIS 



A. Analysis of the S*— wave charmonium decays 

In the previous section, we derived the S*— wave charmonium production rates in 
B decays to Olv"^) for the color-singlet contributions, including one of the color-octet 
contributions at 0{v^). The results depends on three nonperturbative parameters, 
(0|Oi(3^i)|0), (0|Pi(3Si)|0), (0|O8(^5i)|0) or equivalently, |i?v(0)P and Re(R^Vm^) for the 
first two color-singlet matrix elements. Since the third parameter (0|O8(^5'i)|0) is fixed from 
the fit to the ip production at the Tevatron 0, we consider the other two parameters in 
this subsection. In order to determine these two parameters, one has to invoke the lattice 
calculations, some potential models. Or, one can simply determine these parameters from 
the well measured decay rates of J/ip, rj^ which depend on the same parameters. We choose 



6 



the last method to fix two nonperturbative parameters, (0|Oi(^5'i)|0), (0|Pi(^S'i)|0), in this 
work. For this purpose, we hst the decay rates for ip ggg + gg'-f LH{ligh.t hadrons) 
and ^ ^ l+l- ini : 



r(V' ^LH) = 2 Im/i(3^i) {J/i;\0,CSi)\J/^) + 2 Img.CS,) (J/^|Pi(=^^i)| J/^) + 0{v'T), (27) 
r(^ ^ e+e-) = 2 Im/ee(=^5i) {J/^\0,eS,)\J/ij) + 2 ImgU'S,) {J/^Ij\P,CS,)\J/^) + Oiv^T). (28) 

We need to know the short distance coefficients Im/'s to 0{as), and Img's to the leading 
order only in a^, because of the velocity counting rule in the NRQCD [^. From the results 
in the earlier literatures |irT| ||14|, one can extract 



Im/i 

lmgi.^{^Si) 

Imgee^Si) 
Imgqqi^Si) 



lm/3g(J/^ ggg) + lm/i^(J/^ ^ gg-f) 



(vr^-9)(iV,^-4)C^ 
54iV, 



l + (-9.46(2)C. + 4.13(17)C.-hl61(2)„,)^ 



-4 X 4.33 



9)(iV2-4)C, 



2(7r2 - 9)CFQ^a ^ 



4C, 



3A/', 

as 
? — 



3iV, 
-4 X 4.33 

3 

_4 Trg^g^ 
3 3 ' 



54Ar, 
a^(M) 



a.. 



1 + (-9A6Cf + 2.75Ca - 0.77Anf) — 



TT 



2(7r2 - 9)CFQ^a ^ 



af(M), 



13 as 

Ly F 

4 TT , 



(29) 

(30) 
(31) 
(32) 
(33) 
(34) 

(35) 
(36) 
(37) 



In the above expressions, = Ca = Uf = 3 and Cf = 4/3, and Q is the electric charge 
of a heavy quark in the unit of the proton charge. The strong coupling constant as{M) is 
defined in the modified minimal subtraction scheme (MS) for QCD with nj light quarks, 
renormalized at the scale fi = M. The last two are relevant to J/ ijj ^ 'y* ^ hadrons. 
Similar expressions for the rj^ decays are 

r(r/e ^LH)=2 ImfiCSo) {vc\OiC So)\Vc) + 2 lmg^CSo) {vc\PiC So)\Vc) + 0{v'T), (38) 
r(r/c ^ 77) = 2 lmf,^CSo) {Vc\OiCSo)\Vc) + 2 lmg,,CSo) {vc\PiCSo)\Vc) + Oiv'V), (39) 



where 



7CC F 9 
r- 

1 + 



ai{M) 



TT 



--5 



199 

Is 



ISvr^ 
~24 



9 ^ 



(40) 
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Img.CSo) = -l^ (41) 



Im/^^(^5'o) = ttQ 




(42) 



lmg,,{^So) = 7iQ^a\ (43) 



The PDG lists the measured data for these decays as follows |T5| : 



T{il) LH) = (60.72 ± 1.72) keV, (44) 

r(^ ^ e+e") = (5.26 ± 0.37) keV, (45) 

r(r/e ^ LH) = I0.3l:li MeV, (46) 

r{Vc - 77) = 7.0+?;? keV. (47) 

Since the data on rjc decays are not precise enough yet, we use the data on J/ip decays only 
in order to determine two parameters {J / ip\Oi(^ Si)\J / ip) and {J / ip\Pi{^ Si)\J / ip) . Since we 
don't have enough inputs available at this level, we use as{M) to be 0.25 ~ 0.28 instead 
of treating it as a free parameter. Although this choice is not fully systematic from the 
view point of perturbative QCD, our numerical results presented below should not strongly 
depend on the exact value of as{M). From ip LH and ip e~^e~, we determine 

{J/ilj\OiCSi)\J/^) ^ 0.440 (0.490) GeV^ (48) 
{J/ip\Pi{^Si)\J/ilj) ^ 0.025 (0.031) GeV^ (49) 

for as = 0.25 (0.28), respectively. For ijj' , we get 

(7/^'|0i(=^Si)|^') ^ 0.177 (0.198) GeV^ (50) 
(^'|Pi(^'5i)|^/'') ^ 0.008 (0.011) GeV^ (51) 

Note that the radiative corrections to the J/ ijj and ^p' decays are fairly large with or without 
the (0 1 Pi (^5*1)10) term. Radiative corrections increase {J/ip\Oi{^Si)\J/ip) in Eq. (9) by 
~ 80%. The relativistic correction term, {J/ip\Pi{^Si)\J/ip), is about ~ 9% of Eq. (9), 
hence decreases the B ^ J/ip + X rate (Eq. (4)) by ~ 2%. 

Next, let us determine {r]c\Oi{^So)\ric) which is expected to be 

{v,\O^CSo)\vc) = {JMO^CS,)\Jm (1 + 0{v')), (52) 

due to the heavy quark spin symmetry. Since (0|Pi('^S'i)|0) is independent of the total spin 
S, one can determine {ric\Oi{^So)\ric) from one of the decays, Eq. (38) or Eq. (39), and then 
predict the other and compare with the measured rate. Also, the relation (52) should be 
respected in order to be consistent with the NRQCD and the heavy quark spin symmetry. 
If we use (46) as an input, we get (with = 0.25) 

{Vc\OrCSo)\Vc) ^ (0.149;n46) GeV^ (53) 
r{Vc - 77) = (2.8+}:^) keV, (54) 
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the former of which severely violates the heavy quark spin symmetry, (52). For as = 0.28, 
the numbers above change into (^0.1161qo35^ GeV^ and (^l.Sioir) keV. On the other hand, 
if we use (47) as an input with with as = 0.25, we get 

{Vcm'So)\vc)^{0.326t'oZ)G^y', (55) 
T{7], ^ LH) = {23tl) MeV, (56) 

For as = 0.28, the numbers become (OMlto fsl) GeV^ and (32l^) MeV, respectively. Now, 
the relation (48) is better obeyed, although the difference between {J/ ip\Oi{^ Si)\J/ ip) and 
{Vc\Oi{^ So)\ric) is not that small, and the predicted rate for rjc — s>LH is quite large compared 
to the data, (46). One may conclude that the factorization formulation by BBL in terms of 
the NRQCD predicts rather large value of rjc — *>LH, compared to the current experimental 
values, considering the large uncertainties in the measurements. A better determination of 
r(^c — LH) would test the validity of the factorization approach to O(f^). 



B. Results for B decays 

Since all the relevant nonperturbative parameters are in hand now, we are ready to 
estimate the branching ratio for B ^ J/ip + X which includes 0(i;^) corrections and one 
of the O(f^) color-octet contribution. Adding up the change in (OlOf^(^S'i)lO) and the 
color-octet contributions, we get a moderate increase in the branching ratio by factor of 
^ (1.0 + 0.80 - 0.02 + 0.76) = 2.54 of the lowest order prediction (13) for the B ^ J/ijj + X 
case, and (1.0 + 0.80 - 0.02 + 1.14) = 2.92 for the ?/^' case : 

B{B J/tfj + X) = 0.58 %, (57) 
B{B ^x/j' +X) = 0.23 %, (58) 

compared to the data, (0.80±0.08)% and (0.34±0.04±0.03)%. We note that the agreements 
between theoretical estimates and the data get improved, after the radiative corrections and 
the color-octet mechanism have been included. There is a residual uncertainty related with 
the b quark mass Mb. In this work, we have chosen Mb ~ 5.3 GeV, and normalized the 
decay rate to that of the semileptonic B decay in order to reduce the uncertainty from less 
known p. If we use Mb = 4.5 GeV, for example, all the decay rates should be multiplied 
by a factor of 5.3/4.5 ~ 1.18. For B ^ t]c + X , our prediction is tampered by less known 
(0|0^''("'^S'o)|0) as well as (0|Og''(^S'o)|0). Still, we expect that the decay rate increases by a 
factor of ~ 2 or more over the lowest order result, (14). 

Let us finally consider the polarization of J/ip^s produced in B decays. It is convenient 
to define two parameters, ( and a as follows : 

^ TT+L{B^J/^P + Xy ^^^^ 
3C-2 

a = 

The quantity a can be readily measured through the polar angle distribution of dileptons 
in J/if) ^ in the rest frame of J/ip : 
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(61) 



where 9 is the angle between the flight direction of a lepton in the rest frame of J/V' and the 
flight direction of J/ip in the rest frame of the initial B meson. For the unpolarized J/'ip, 
we would have C = 2/3 (a = 0) corresponding to the flat cos6' distribution of dileptons. 
Assuming the factorization for B ^ J + X in the color-singlet component in the effective 



Hamiltonian (2), one flnds that J/i/j^s produced in B decays are polarized with [|T6 



a = -—I -I ^ -0.49 (-0.36) (62) 

or C = 0.41 (0.49) for = 5.3 (4.5) GeV, which nicely compares with the CLEO 
measurement ae^p = —0.44 ± 0.16. One may wonder if the color-octet mechanism 
considered in this work can change the polarization of J ftp substantially. However, the 
structure of the amplitude for b ^ J/ip + X due to the color-octet ^Si state is the same as 
that due to the color-singlet mechanism (including both the lowest and the next-to-leading 
order terms in f ^). Namely, the amplitude for b ^ J/ip + X is proportional to 

Mcxe^sY{l-l5)b. (63) 

Thus, the prediction for a remains the same even in the presence of the {cc)s(^Si) color-octet 
contribution. This is in sharp contrast with the case of ip production through the gluon 
fragmentation into the color-octet cc state which in turn evolves into ip . In the case of 
gluon fragmentation {g {cc){^Si)s J/ip + X), the initial gluon (with ^ (2Mc)^ and 
g° >> 2Mc) is almost on shell, being almost transverse up to /qq- Thus, the polarization 
of the color-octet cc is also almost transverse. Because of the spin symmetry of heavy quark 
system, the polarization of the daughter J /if) is the same as the parental color-octet cc state. 



IV. CONCLUSIONS 

In this work, we considered the relativistic corrections and a color-octet contribution to 
S'— wave charmonium productions in B decays. Our results, (57) and (58), give moderate 
increase to the previous analyses based on the color-singlet mechanism in the nonrelativistic 
limit. Compared to the previous analyses of the lowest order in v"^ and a^, we get an ~ 80% 
increase in [J / %Ij\Oi(^ Si)\J / ip) from the radiative corrections to J/ip — > e'''e~ and J/ip — *• 
LH, ~ 2% decrease from the relativistic correction through the {J / %Ij\Pi{^ Si)\J / ip) term, 
and ~ 76% increase in the decay rate from the color-octet contribution, (0|O8^^(^5'i)|0). 
Thus, the color-octet mechanism, which has been proposed as a possible solution to the ip 
puzzle at the Tevatron, could give an enhancement of the B ^ J/ip + X decay rate by a 
moderate amount. 

It should be kept in mind that what we considered in this work is only one of the color- 
octet operators that may contribute to B ^ J/tp + X. We have chosen only {0\Oi^^{^Si)\0), 
since we know the numerical value of this matrix element from the work of Braaten and 
Fleming ^ . This matrix element is rather special in the sense that it is the only color-octet 
operator which is relevant to g (cc)8(^5'i) ^ J/ip + X in the leading order in v"^ and a^. 
However, for the B decays, other color-octet operators can contribute as well ; e.g.. 
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r{B^{cc)8{ Si)+X ^7]c + X) = {C+ + C-) + ^0, (65) 

and similar expressions for the contributions of {0\Og{^Pj)\0) : 

TiB ^ (cc)8(^P0 +X^J/i; + X)= (0|Qb^24^i)|0) ^ ^_)2 (^1 + ) To, (66) 

TiB iccU'P,) +X^r^, + X)= (C+ + C.r (l + ^ ) To, (67) 

In order to estimate the effects of these color-octet matrix elements (0|O8'''^( ^5*0)10) and 
(0|O8'''^('^Pi)|0), we need to consider other processes as well, such as 7 + p J/ip + X and 
e~^e~ — s> 7* — > J/V^ + X For example, the height of the elastic peak for photoproduction 




of J/ip depends on these color-octet matrix elements 0. Complete analysis of color-octet 



contributions to the 5*— wave charmonia in 7p collisions is called for as well ||19[. Once 
these new color-octet matrix elements are determined from other processes, our results in 
this work will provide an independent test of the hypothesis of color-octet mechanism as a 
possible solution to ip anomaly at the Tevatron. 

We have also analyzed the leptonic and the inclusive hadronic decays of J/ijj and ijj to 
O(f^) in the framework of the BBL's factorization scheme, and did extract two nonpertur- 
bative parameters, {H\Oi{^Si)\H) and {H\Pi(^Si)\H) with H = J/ip and ip . These are 
important inputs in many other theoretical calculations of the S*— wave charmonia produc- 
tions in various high energy processes and their subsequent decays. As a by-product, we 
have found that the inclusive hadronic decay rate for 77^ may be larger than the current 
PDG value by factor of ~ 2, if the BBL's factorization formulae to O(w^) works with the 
charmonium system. The better measurements of T{ric LH) would test our predictions 
based on the factorization approach for the heavy quarkonium decays in the framework of 
NRQCD to 0{v^). 
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